Intensity distribution of scalar waves propagating in random media 
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Transmission of a scalar field through a random medium, represented by a system of randomly 
distributed dielectric cylinders, is calculated numerically. The system is mapped to the problem of 
electronic transport in disordered two-dimensional systems. Universality of the statistical distribu- 
tion of transmission parameters is analyzed in the metallic and localized regimes. In the metallic 
regime, the universality of transmission statistics in all transparent channels is observed. In the 
band gaps, we distinguish a disorder induced (Anderson) localization from tunneling through the 
system, due to a gap in the density of states. We also show that absorption causes a rapid decrease 
of the mean conductance, but, contrary to the case of the localized regime, the conductance is 
self-averaged with a Gaussian distribution. 
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PACS numbers: 41.20.Jb,72.10.-d, 73.23.-b 



I. INTRODUCTION 

Transport of classical waves in random media is a chal- 
lenging problem, attracting increasing interest of theo- 
retical and experimental physicists because it offers the 
possibility of studying Anderson localization.^ Since in- 
teractions do not play any role in classical wave scat- 
tering, these systems might be more convenient for ex- 
perimentally verifying the scaling theory of localization^ 
than quantum electronic systems, where the influence of 
the mutual interaction of electrons upon transport has 
not been yet clarified>^ Two main issues of localization 
theory, namely the presence or absence of the metallic 
state in the two-dimensional (2D) systems, and the va- 
lidity of the single parameter scaling (SPS)^'"*'^ might 
be more readily resolved experimentally for the classical 
wave problem than for the electronic one. Recent exper- 
imental results for the transmission of electromagnetic 
waves indeed confirmed that transmission is universal in 
the diffusive regime,^ and presented strong indications 
for disordered induced Anderson localization. ^'^'^ 

In this paper, we analyze numerically the transmis- 
sion of scalar classical waves through a two dimensional 
(2D) system of randomly distributed dielectric cylinders. 
Following Ref. we map the problem into the 2D An- 
derson model with random binary potential. Statistical 
properties of wave transmission are then analyzed using 
the transfer matrix method;— We calculate the conduc- 
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In Eq. Tab = l^abP, where tab is the transmission am- 
plitude from channel a to channel b. a, 6 = 1, 2 . . . Nop, 
where N^p is the number of open channels. We first de- 
termine the band structure of the original classical wave 
problem. Then, we analyze the statistical properties of 
the transmission in bands and gaps. 

In bands, where g > I, we observed diffusive trans- 
port. Statistical properties of the transmission are in 



good agreement with theoretical predictions of the ran- 
dom matrix theorjiMii^ and the DMPK equation>iSiii 
The distribution of the conductance is Gaussian with a 
universal dimension-dependent variance li^siSiiSiSSiSi Uni- 
versal properties were predicted not only for the conduc- 
tance g, but also for the normalized parameters Sab = 
Tab/{Tab)^^^ and for the normalized transmission in a 
given transport channel 
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The universal probability distribution 
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has been derived analyticallyi^^iS^ From Eq. Q , the sec- 
ond cumulant is obtained as 
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Universality of the statistical properties of parameters, 
Sa, was confirmed experimentally^ up to rather small val- 
ues of the conductance (((/)« 2 — 3). 

In gaps, the mean of the logarithm (In g) decreases lin- 
early with the system size. Here, we distinguish between 
two different regimes, one with non-zero density of states, 
and the other called tunneling regime, characteristic for 
the frequency region without eigenstates. The first de- 
scribes Anderson localization, characterized by the Gaus- 
sian distribution of (In g) with variance var (In oc —L/S, 
(^ is the localization length) in agreement with local- 
ization theory. The second regime appears in the gaps, 
where the density of states is very small. Then, the trans- 
mission is determined by tunneling through the sample. 
Although (In g) cx — L, the distribution of (In g) is not 
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FIG. 1: Typical sample configuration. The size of the system 
is 128 X 128, the filling factor x = 0.1. 182 rectangular rods 
with square cross-section of 3 x 3 are randomly distributed 
throughout the sample. Scalar wave propagates from left to 
the right. 



Gaussian but given by the statistics of the energy spec- 
tra. 

The form of the probabihty distribution p(ln g) enables 
us also to distinguish between localization and absorp- 
tion. We find that absorption also gives a decrease of 
(In g) cx — L; however, in contrast to localization, the 
conductance remains self-averaged. Our data agree with 
the theoretical predictiouc^Si as well as with the qualita- 
tive criterion^ for localization 2/(3varSa) < 1. 

The paper is organized as follows: In Sect. ^ we 
introduce the model and system parameters. In Sect. IIIII 
we present the mapping of the classical wave problem into 
the problem of the transmission of electrons in disordered 
systems^ Results of numerical simulations are presented 
in Sect. IIVI and Ivl Conclusions are given in Sect. IVII 



II. MODEL 

We study a two-dimensional system consisting of a ran- 
dom array of dielectric cylinders. To make the numerical 
simulations easier, we consider rectangular rods instead 
of cylinders. A typical sample is shown in Fig. ^ No 
contact or overlap of neighboring rods is allowed. The 
dielectric permittivity of rods is £2 and of the embedding 
medium is ei. Two semi-infinite leads with permittivity 
eicad are attached to the sample. The concentration of the 
dielectric rods is given by the filling factor x of the rod's 
material. In this work we consider square cross-section 
of dielectric 3x3 rods, measured in dimensionless units. 
In the same units, the system size L varies from 32 up to 
256. 



III. MAPPING TO THE ELECTRONIC SYSTEM 

We use the formal equivalencsiS. of the scalar wave 
equation, 

\/^u + e{f)uj^u^O, (5) 

and the Schrodinger equation for electrons, 

V^* + [E- V{f)] ^' = 0, (6) 

(we set 2mc^/h^ — 1 and = 1). For a given space 
dependence, e(r)j and given frequency, w, one can find 
a potential V{f) and energy E such that the solutions u 
and 5* of Eqs. (O and (0 are identical. Formal equiva- 
lence requires that the identity 

coheir) = E- V^r) (7) 

must be fulfilled for all f. 

Formula |7J| does not mean that the two models de- 
scribed by Eqs. Q and ® are equivalent, because po- 
tential V{r) depends both on the frequency, w, and the 
energy, E. Equation Q only means that for a given 
scalar wave model defined by e(r), we can find for each 
frequency, uj, an electronic model with potential V{r) 
and energy E, that the solutions of both models are the 
same. 

Both the energy and the potential of the electronic 
model are fixed by formula Q) . Changing the frequency, 
UJ, we obtain another electronic model, since the poten- 
tial V{r) changes. This means that for a given spatial 
distribution of e{r), two different frequencies, uj, define 
two different electronic models. 

To be more specific, we consider a model of randomly 
distributed rectangular rods discussed in Sect. ^1 This 
model can be mapped into the electronic model with a 
random binary potential. For a given sample, the spatial 
distribution of the potential is identical with the distri- 
bution of the permittivity. The energy, E, as well as two 
values of the potential, Vi and V2, is determined by Eq. 
([7|) with the following relations 

w'eiead = E (8) 
Lo^ei ^ E-Vi (9) 
u;h2 = E~V2. (10) 

From Eqs. © and H1U|) we easily obtain 

^'e, = (11) 

where 

/^=- (12) 
ei 

and 

S = V2-Vi, (13) 
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The energy, E, is given as 
E = Vi 
and Vi is determined by 



1' 



Vi = uj^{ei - eicad)- 



(14) 



(15) 



Note that for ejoad = = 1 (vacuum in leads), Eqs. 
and ((ini can be solved only when both Vi, V2 > 
assuming that ei, £2 > 1. 

For simplicity, we consider in this paper the special 
case, 



(16) 



The method is, of course, applicable to any set (eicad, ei, 
62), including eicad > ei- 

Using Eqs. p5|l and H16|l. we have Vi = 0. Comparison 
of Eq. ini) and Eq. lO gives lu'^ = E. Finally, V2 = -6. 
We also fix the ratio of two permittivities fi = €2/^1, to 
the value of /i = 11. 

In numerical simulations, we used the discretized ver- 
sion of the Schrodinger equation, Eq. 
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(17) 

Hard wall boundary conditions were used. The transfer 
matrix method^i was used to calculate numerically all 
the parameters, Sa, and the conductance. The number 
of open channels, A'op < L, which enters in Eq. ^ is 
given by the number of propagating solutions (A:„ real) 
of the dispersion relation 



2 cos kn ^ 4: — E — 2 cos ■ 



1 



1,2,. ..,L. (18) 



Each frequency, uj, defines a corresponding tight-binding 
Hamiltonian, for which the transmission is calculated us- 
ing standard numerical procedures known for the elec- 
tronic localization problems,-?^- Statistical ensembles of 
-^stat — 10^ samples were considered, which assure that 
sufficient accuracy of the transmission parameters of in- 
terest were obtained. 



IV. CONDUCTANCE 

Figures |2| and O show the oj dependence of the mean 
conductance for the present model. Equivalently, the nu- 
merical results of Figs. |21 and |21 might be interpreted as 
the (5-dependence of the conductance of the 2D electronic 
system. In the last case, however, one must keep in mind 
that in the Anderson model, the value of 5 determines 
the random potential. Data presented in Figs. [21 and 
Oldo not directly correspond to the electronic density of 
states in the disordered electronic system. Different val- 
ues of i5 correspond to different models. Note also that 
the energy, E, Eq. (|14|l is also a function of S. 

For completeness, we show in Fig. 0] the density of 
states, p{E), of the electronic system for four values of 6. 
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FIG. 2: Mean conductance (g) (top) and (ing) (bottom) as a 
function of the frequency for /i = £2/61 = 11 and x = 0.2. Fre- 
quency bands and gaps are clearly visible. In bands around 
LJ^ — 0.14 and 0.31, we find transport statistics typical for 
the metallic regime. Two different transport regimes were 
observed in gaps: In the first gap, a;^ ~ 0.08, the density of 
states is very small and transport is due to tunneling through 
the sample. In the second gap, lo'^ ~ 0.21, we observed disor- 
der induced Anderson localization. 
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FIG. 3: The same as in Fig. El but for x = 0.1. 

A. Diffusive regime {{g) > 1) 

Three frequency bands are visible in Figs. 13 and O 
where we expect the metallic behavior. Mean conduc- 
tance {g) > 1 slightly increases with the system size and 
the value of var g is very close to the universal conduc- 
tance fluctuationiiSiSSiSi The distribution of the conduc- 
tance is Gaussian (see, for example. Fig. Ill|l . Although 
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FIG. 4: Density of states of the electronic system with x — 0.2 
calculated for four values of 5. The size of the system is 48 x 48 
and an average over 10 ensembles was calculated. Dashed 
lines indicate the energy, E{S). For S — 0.8, the inset shows 
the position of eigenenergies for 18 different samples. 




FIG. 5: Left panels: Probability distribution p{sa) for some 
open channels compared with the theoretical prediction of 
Eq. (EJ (dashed line), uj^ = 0.14 (the center of the first 
band). Concentration of cylinders is a; = 0.1 (top) and x — 0.2 
(bottom). Right panels: The second cumulant ^{g)varsa as a 
function of index a (top) and mean values (Ta) (bottom). The 
size of the system is L = 192, electron energy is E = 0.14, 
the number of open channels, A*'op — 23. Statistical ensemble 
of A'stat ~ 10000 samples was considered. 
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FIG. 6: The same as in Fig. 1^1 but for lj^ = 0.31. iVop = 
34. The agreement with theory is not as good as in Fig. |3 
especially for a: = 0.1 because of the small system size. 
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FIG. 7: The second cumulant | {g)va,TSa as a function of index 
a for the center of the second band u)'^ — 0.31 and a; = 0.1. 2D 
systems of size L x L with L — 125, 192 and 256 were consid- 
ered. These results are compared with quasi-one dimensional 
systems of size 128 x with = 128, 256, 320, 448 and 512 
with mean conductance (g) = 7.28, 4.29, 3.53, 2.57 and 2.25, 
respectively. 



these properties are finite size effects (no metallic state 
exists in 2D in the hmit L — > oo), the numerical data 
enable us to check the theoretical prediction about the 
transmission statistics. 

In Fig. |31 we show the statistics of the parameters, 
Sa, for the frequency in the center of the first {S « 1.4) 
band. The number of open channels is 23 for the size of 
the system 192 x 192. Results confirm that the distribu- 
tion, P{sa), is universal and does not depend on a . The 
second cumulant, varsa, is close to its theoretical value, 



2/3(5)- 

The same analysis was completed for the second band 
(Fig. EJ. Here, the agreement with theory is not as good 
as in the previous case, presented in Fig. (Sjespecially for a 
smaller concentration of rods. Although the distribution, 
P(sa), does not depend on a, it differs considerably from 
theoretical predictions. We interpret this discrepancy as 
a finite size effect. Indeed, the mean conductance (g) 
increases with the system size (Fig. O and inset of Fig. 
|S} , which indicates that we have not reached the diffusive 
regime yet. 
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FIG. 8: The variance of the conductance, var g = {g^) — {g)'^, 
as a function of the mean conductance for a; = 0.1 and x = 0.2, 
and various system sizes. Open symbols: x — 0.1, full sym- 
bols: X — 0.2. For small values of {g), the numerical results 
scale to a universal curve. For large (g) they converge to the 
universal value of var q = 0.1855^^ The inset shows the 
mean conductance {g) for S — 1.4 (circles) and 5 — 3.1 (trian- 
gles). For S = 3.1, {g) still increases with L, which indicates 
that L is not large enough for transport to be diffusive. 
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FIG. 10: The probability distribution p{ln g) for u;'^ = 0.08 
and X = 0.2. The distribution shows that transport is due 
to tunneling through a system with an energy gap. Note 
that the width of the distribution, var Ing -C (Ing), is much 
smaller than that of the disorder induced insulator. Data for 
the mean value and variance of In g are given in the legend. 
Nst^t = 7000 for L < 256 and is > lO" for smaller L. 



B. Localization ({g) <^ 1) 



L <ln g> var In g 




FIG. 9: The probability distribution P{lng) for uj'^ = 0.21. 
The distribution is close to the log normal with var Ing « 
— (lug), typical behavior of Anderson localization. The num- 
ber of random configurations A'stat is 10^ for L < 96 and 10* 
for L = 128. The legend presents the data for (Ing) and var 
Ing. 



Note also that we are studying 2D samples, while the 
theory is formulated for quasi one dimensional systems. 
Therefore, we expect that the agreement with theory 
should be better if the length of the system increases. 
This is confirmed by numerical results presented in Fig. 
13 which shows the second cumulant ^{g)va.T Sa for vari- 
ous 2D and quasi one dimensional system. 

Figure 121 shows the variance, var g, of the conductance 
vs the mean conductance (g). We expect that var g 
should saturate to the universal value of var g 0.1855 
[18.2O.J for 5 > 1. 



In the regions between the pass bands, the conductance 
decreases exponentially with the system size. It is much 
more pronounced for larger concentrations of cylinders 
(cc = 0.2). Here, two different regimes were observed. 

In the upper gap, {6 « 0.21) the observed statisti- 
cal properties of the conductance are in agreement with 
the theoretical expectations for the localized regime: 
The distribution of \ng is Gaussian (Fig. O with var 
ln(7 « —(In (7) oc 2L/^, which is characteristic for the 
disorder induced localization. The parameter ^ is the 
localization lengthi^ 

In the lower gap {S « 0.08), the conductance decreases 
rapidly as the size of the system increases. The proba- 
bility distribution of In g is not Gaussian, as can be seen 
in Fig. ^1 Instead, it decreases exponentially for larger 
values of conductance as 

p{\iig) oc expconst[(ln5) — lug]. (19) 

In contrast to Anderson localization, no samples with 
conductance close to 1 were found. This indicates that 
transport is possible only by tunneling through isolated 
eigenstates. Also, the variance, var In g, is much smaller 
than that of the Anderson insulator, 

var Ing < -(Ing). (20) 

The last property seems to be in agreement with previous 
work of Deych et alm^ who argued that single parameter 
scaling does not work in the energy intervals, where the 
density of states is so small that another characteristic 
length Is sin^ ^ p{E) exceeds the localization length. 
Indeed, we found that the density of states is close to 
zero in the neighborhood oi E = 0.08 (Fig. 0)), so that 
the average distance between isolated eigenstates is larger 
than the localization length (^ « 8 was estimated form 
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FIG. 11: The probability distribution of the conductance g 
for X = 0.2 and lo"^ = 0.14 for the case without absorption and 
with small absorption in the dielectric cylinders (L — 128). 
Absorption is due to the small imaginary part of the permit- 
tivity £2 of cylinders. The mean conductance decreases, due 
to absorption (see legend) and the conductance distribution 
is Gaussian. The inset shows the L dependence of \n{g) and 
(Ing) for the stronger absorption (Im V2 = 0.05). The numer- 
ical results confirm that the conductance is self-averaged. 



the L - dependence (Ing)). For instance, in the inter- 
val 0.05 < E < 0.11, we found an average of only one 
eigenstate, when L = 48 so that Ig > 48. Although Ig 
is expected to decrease when L increases, we were not 
able to reach Anderson localizatfon, even for the largest 
system studied, L = 256. 



V. ABSORPTION 

Absorption reduces the transmission of the EM waves 
in a similar way as localization. Mean conductance de- 
creases exponentially with the system length;--. To dis- 
tinguish between Anderson localization and absorption 
effects, we need to understand the statistical properties 
of the transmission. In Ref. the simple criterion for 
localization was derived, based on the value of the pa- 
rameter 



g' = 2/(3varsQ). 



(21) 



It was argued that localization appears if g' < 1. 

To study the effects of absorption, we add a small 
imaginary part to the permittivity of cylinders (more 
precisely, V2 in Eq. IjlOfl becomes complex in our sim- 
ulations). First, we analyze how absorption changes the 
transmission properties of the metallic system. We use 
5 = 1.4 and x = 0.2 (Fig. [^J. As expected, the mean 
conductance decreases when the system size increases, 
(In g) ~ —L similarly as for localized waves. In contrast 
to the localized regime, the conductance is self-averaged 
in this case, (In g) = ln(g). The conductance distribution 
is still Gaussian. As shown in Fig. ^2 the width of the 
distribution of the normalized conductance depends only 
weakly on the absorption strength. This is in agreement 
with analytical results of Brouweri^ 




FIG. 12: The probability distribution of the logarithm of the 
conductance in various regimes. L: regime of Anderson lo- 
calization with Gaussian distribution and var In y « — (In g) . 
LA: localized regime with absorption (Im V2 = 0.01). Com- 
pared with L, we see that the part of the distribution with 
relatively large conductance is missing. T: tunneling regime, 
in which varln(/ <^ —(Ing). TA: tunneling with absorption 
(Im V2 = 0.05). Absorption does not influence p{lng). MA: 
metallic regime with absorption (Im V2 ~ 0.065). Here, con- 
ductance is self-averaged. Parameters of systems were chosen 
such that (Ing) is approximately the same for all systems 
(see legend). Inset shows var Sa in all the above regimes. As 
predicted in Ref. |^ var Sa < 1 in the metallic regime with 
absorption, but is ^ 1 in the localized regime, both with and 
without absorption. The same holds for the tunneling regime 
(data for L, LA and TA regimes are almost indistinguish- 
able). Thus, the criterion g' < I cannot distinguish between 
localization and tunneling. 



Figure El compares the statistics of the logarithm 
of the conductance of five different transport regimes 
~ localization and tunneling with and without absorp- 
tion, and diffusive (metallic) regime with absorption. 
The numerical results confirm that absorption does not 
change the statistical properties of the conductance, 
given mostly by the interference of electrons or classical 
waves due to disorder. In the regime of Anderson local- 
ization, the probability to find relatively large values of 
g is reduced due to absorption, while another part of the 
distribution, where Ing (In g) , is almost unaffected by 
the presence of absorption. As a result, pilng) is not 
Gaussian anymore. In the tunneling regime, absorption 
only reduces the magnitude of var Sa as can be seen in 
the inset of Fig. [T3 



VI. CONCLUSIONS 

We presented a numerical analysis of the transmission 
of the scalar classical wave through a disordered two 
dimensional system. By mapping this system to a 2D 
electronic disordered problem, we found frequency inter- 
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vals with high transmission. For these frequencies, we 
obtained the statistical distribution of the transmission 
parameters predicted recently by the theory. We con- 
firmed the universality of the conductance fluctuations 
and of the distribution of the parameters Sa in the metal- 
lic regime. The universality survives for rather small val- 
ues of conductance, (g) = 1 — 3. Our numerical results 
confirm that the theory, developed for the quasi one di- 
mensional systems, can be successfully applied to 2D sys- 
tems, too. Our results are also consistent with previous 
experiments]^ 

As there is no metallic regime in two dimensional sys- 
tems, the above described metallic behavior is just an 
effect of the finite size of our sample. By increasing the 
system size, conductance would decrease and finally, in 
the limit of L ^ localization length, the wave becomes lo- 
calized. The localized regime was also observed in a gap, 
where the density of states is smaller than in the bands, 
but still non-zero. In this frequency interval, broad dis- 
tribution of the logarithm of the conductance, typical for 
the Anderson localization in electronic systems, is ob- 
served. 

Anderson localization should be distinguish from the 



tunneling regime, which we found in the frequency gap, 
where the density of states is close to zero. Here, the 
statistics of conductance is determined by the statisti- 
cal properties of the isolated frequencies inside the gap. 
The distribution p{\ng) differs considerably from Gaus- 
sian. These results are in agreement with the theoretical 
analysis of Deych et alm^ 

Finally, we analyzed the effects of absorption. We 
found, in agreement with theoretical and experimental 
works, that absorption does not change the statistical 
properties of the transmission. Since the statistical prop- 
erties of the parameters, Sa, are insensitive to the pres- 
ence of absorption both in the localized and in the metal- 
lic regimes, typical values of var Sa enable us to decide 
whether the exponential decrease of the conductance is 
due to localization or to absorption. However, statistics 
of Sa can not distinguish between Anderson localization 
and tunneling. Fortunately, these two regimes are distin- 
guishable from the form of the probability distribution of 
the total transmission g. 
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